Demonstration of a quantum logic gate in a cryogenic surface-electrode ion trap 
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We demonstrate quantum control techniques for a single trapped ion in a cryogenic, surface-electrode 
trap. A narrow optical transition of Sr^ along with the ground and first excited motional states of 
the harmonic trapping potential form a two-qubit system. The optical qubit transition is susceptible 
to magnetic field fluctuations, which we stabilize with a simple and compact method using super- 
conducting rings. Decoherence of the motional qubit is suppressed by the cryogenic environment. 
AC Stark shift correction is accomplished by controlling the laser phase in the pulse sequencer, elim- 
inating the need for an additional laser. Quantum process tomography is implemented on atomic 
and motional states by use of conditional pulse sequences. With these techniques, we demonstrate 
a Cirac-ZoUer controUed-NOT gate in a single ion with a mean fidelity of 91(1)%. 



I. INTRODUCTION 

Trapped ions are promising candidates for realizing 
large-scale quantum computation P, Significant 
progress has been made in demonstrating the funda- 
mental ingredients of a quantum processor, with much 
progress in gate fidelities ^] and multi-ion entanglement 
[J, [^]. In recent years there has been increasing inter- 
est in microfabricated surface-electrode traps, owing to 
their inherent scalability [1, However, quantum gates 
have yet to be demonstrated in such systems. An issue 
with miniaturization of traps is that anomalous heat- 
ing of the ion's motional state scales unfavorably with 
trap size potentially limiting gate fidelity in traps 
of suitable dimensions for scalability [9||. Recently, it 
has been shown that by coofing to cryogenic tempera- 
tures, the heating rate can be reduced by several orders 
of magnitude from room-temperature values [To| . thus 
providing one potential solution to this problem. In this 
work, we demonstrate a quantum gate in a microfab- 
ricated surface-electrode ion trap that is operated in a 
cryogenic environment, and present some control tech- 
niques developed for this experiment. 

We implement a Cirac-Zoller controUed-NOT (CNOT) 
gate using qubits represented by the atomic and motional 
states of a single ion. The S ^ D optical transition in 
^^Sr+ is used as one of the qubits. The motional ground 
state and first excited state of the ion in the harmonic 
trap potential form the second qubit. The optical tran- 
sition has the advantage of a long lifetime while requir- 
ing only a single laser (unlike hyperfine qubits), but the 
qubit is first-order Zeeman sensitive, which makes it sus- 
ceptible to magnetic field noise. Taking advantage of the 
cryogenic environment, we stabilize the magnetic field us- 
ing a pair of superconducting rings Since the Sr+ 
ion qubit is not an ideal two-level system, the coupling 
between the sideband and carrier transitions causes level 
shifts known as the ac Stark shift, which must be cor- 



rected. In previous work, this has been accomplished 
with an additional laser field with the opposite detuning 
to cancel out the shift [12]. Here, to reduce the experi- 
mental complexity of the additional acousto-optical mod- 
ulators (AOMs) and optics required. Stark shift correc- 
tions are implemented in the experiment control scheme 
by shifting reference frames as is done in NMR [l^ . For 
readout, the qubit encoded in the motional state of the 
ion normally cannot be measured directly, but condi- 
tional pulse sequences allow full state tomography of the 
qubit system. 

The control techniques developed here may be applica- 
ble to use of a single ion to probe and manipulate other 
systems, even though they focus on a single ion and do 
not necessarily imply scalability. Some such systems in- 
clude the coupling of ions to superconducting qubits [13] , 
micromechanical cantilevers [15|, cavities [IQ] . and wires 
[l7| . In many of these experiments, maximizing the cou- 
pling requires proximity of the ion to a surface, and co- 
herence of the motional state is also desired. 

This paper is organized as follows. The experimental 
setup, including the magnetic field stabilization scheme, 
is described in Sec. [Hi Section Hill briefly discusses mo- 
tional state decoherence and shows that such decoherence 
has an insignificant effect on the gate performance in our 
system. Section IIVI presents a theoretical model of the 
Stark shift correction and experimental implementation 
of the method. Section |V] describes the state prepara- 
tion and measurement sequences that allowed us to im- 
plement quantum process tomography on the single-ion 
system. Section lVTl describes the realization of the CNOT 
gate, along with a discussion of gate performance and er- 
ror sources. 



II. EXPERIMENTAL SETUP 



A. Cryogenic microfabricated trap 



1sxwang@mit .edu] 



The microfabricated trap is a five-rod surface-electrode 
design identical in geometry to that described in Ref. 
[13 . The trap is made of niobium, and the fabrica- 
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tion process is similar to prior methods [lO| employed 
for gold traps, and is described briefly here. A 440-nm 
Nb layer is grown on a sapphire substrate by sputtering. 
The sheet resistance is 0.3 fi/sq at 295 K, and the super- 
conducting transition is at Tc = 9.f 5 K. Trap electrodes 
are patterned using NR9-3000 photoresist and etched us- 
ing reactive ion etching with CF4-I-02. The trap center 
is 100 /^m above the surface. For the experiment de- 
scribed here, the axial and two radial trap frequencies are 
2-K X {1.32,2.4,2.7} MHz, respectively Although a su- 
perconducting trap was used for the work described here, 
the effects of the superconducting material on trapping 
behavior will be described elsewhere jl^ . 

The trap is cooled and operated in a 4 K bath cryostat 
described in Ref. [ii]. Typical ion lifetime is on the order 
of several hours, limited only by the liquid helium hold 
time. Loading is done via photoionization of a thermal 
vapor. 
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FIG. 1: (a) **Sr+ level diagram. The 422- and 1091-nm 
transitions are used for Doppler cooling and detection. The 
673.837-nm transition couples the qubit levels, (b) Details 
of the qubit states with Zeeman levels explicitly drawn. 
The "pump" transition is used to pump the ion out of the 
Si/2ijn =1/2) state during initialization. 



B. Sr qubit and laser system 

An atomic ion confined in a harmonic trapping po- 
tential can encode two qubits, one in its optical atomic 
transition and one in its lowest motional states. The 
^^Sr"*" ion has a narrow optical transition, 6*1/2 ^ ^5/2 
with a linewidth of 0.4 Hz. The m = -1/2 o m = -5/2 
levels are used for the atomic qubit transition. This tran- 
sition is chosen for convenience, as along with the 
(to — —3/2) level it forms a closed three-level system for 
sideband cooling [23|. The degeneracy of the multiple 
Zeeman levels is lifted by applying a constant field of 4 G 
with external coils. To address this transition at 674 nm, 
a diode laser is grating-stabilized and locked to an exter- 
nal cavity via optical feedback It is further stabi- 
lized by locking to a high-finesse cavity made of ultralow 
expansion glass as in Ref. [23| . The frequency noise, indi- 
cated by the Pound-Drever-Hall error signal as measured 
with a spectrum analyzer, is 0.3 Hz for noise components 
above 1 kHz. Below 1 kHz, acoustic noise broadens the 
laser linewidth to ^300 Hz, an estimate based on Ramsey 
spectroscopy measurements on the carrier S — D transi- 
tion assuming that the Ramsey contrast decay is caused 
primarily by the laser linewidth. This laser beam prop- 
agates along the axial direction of the trap, so we ig- 
nore the radial modes of motion in sideband cooling and 
quantum operations. Doppler cooling is performed on 
the Si/2 ^ Pi/2 transition with a 422 nm diode laser. 
Two IR diode lasers, at 1092 and 1033 nm, repump the 
ion from the -D5/2 and D3/2 states. For all measurements, 
the ion is initialized to the Si/2{m = —1/2) state and the 
motional ground state via a sequence of Doppler cooling, 
sideband cooling, and optical pumping. Figure [T] shows 
the relevant levels of Sr+ for the experiment. 

A pulse sequencer [2^ consisting of a field- 
programmable gate array (OpalKelly XEM3010-1000) 
and direct digital synthesis boards controls the phases, 
amplitudes, and lengths of the laser pulses. Switching of 



the beam and setting of the desired frequency and phase 
shift are accomplished using AO Ms on the 674-, 422-, 
and 1033-nm lasers. Phase- coherent switching is imple- 
mented by computing the expected phase at time t, refer- 
enced to a fixed point in the past, for a given frequency 
/ using (j)o{t) = ft{mod 2tt). Then, after a frequency 
switch at time T, the absolute phase of the waveform 
is adjusted to equal ipoiT) -I- 0, where </> is any desired 
phase. This process allows for frequency switching while 
maintaining phase information throughout any arbitrary 
pulse sequence. 



C. Magnetic field stabilization 

When the optical qubit is encoded in a pair of levels 
that are first-order sensitive to magnetic fields, field fluc- 
tuations on the time scale of gate operations will decrease 
gate fidelity. One way of passively stabilizing the field is 
by use of a /i-metal shield, which is expensive and incon- 
venient for optical access, and also mainly effective for 
low-frequency noise. Active stabilization of the magnetic 
field using a flux gate sensor and coils has been imple- 
mented in another experiment [IJl , at the cost of higher 
complexity. 

Superconducting solenoids have been employed for 
passively stabilizing ambient magnetic field fluctuations 
in NMR experiments, with field suppression by a factor 
of 156 ^J4|. A similar method for ion traps which would 
permit good optical access is desired. In the NMR imple- 
mentation, the field needs to be stabilized over a region 
1 cm in length, whereas in an ion trap the region of inter- 
est is much smaller. Our method uses the same principle 
of superconductive shielding, but the small region and 
requirement for optical access suggest a more compact 
approach. 

We stabilize the magnetic field by employing the per- 
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sistent current in two superconducting rings, placed 
closely adjacent to the ion trap chip. This is a very com- 
pact and experimentally convenient arrangement, with 
high passive field stability and little barrier to optical ac- 
cess. Below the trap is a 1 x 1 cm^ square Nb plate with 
a 1.5 mm diameter hole, located 0.5 mm below the trap 
center. Above the trap is a 50 cm^ square plate with 
an 11-mm-diameter hole, located 7 mm above the trap 
center [Figl2ja)]. Both rings are 0.5 mm thick. This ge- 
ometry was chosen to optimize the field suppression at 
the trap location using the method to calculate rnagnetic 
fields in superconducting rings described in Ref.[25|. 

With a single trapped ion we measured the field sup- 
pression by applying a constant field with external coils, 
cooling the trap and Nb rings to below Tc, and reduc- 
ing the field while measuring the 5* O D transition fre- 
quency. The magnetic field is calculated from the Zeeman 
splitting between the m = —1/2 0771 = —5/2 transition 
and the m = -1-1/2 ^ m — —3/2 transitions. A 50- 
fold reduction in field sensitivity was observed (Figl2]), in 
agreement with the numerical calculation. To determine 
the effectiveness of the noise suppression on coherence 
of the atomic qubit, we measured the decay of Ramsey 
fringes as a function of the separation of the Ramsey 7r/2 
rotations on the carrier S ^ D transition. Such a mea- 
surement also includes effects caused by laser linewidth 
and the drift in laser-ion distance. We found that reduc- 
ing the magnetic field noise by a factor of 50 did not im- 
prove the coherence time by more than a factor of 2, from 
T2 ~ 350 ^s to ~ 660 fis. This suggests that magnetic 
field noise is no longer a dominant source of decoherence 
when compared to laser linewidth. Although this mea- 
surement was done under dc and the dominant source of 
magnetic field fiuctuations is frequencies near 60 Hz and 
its harmonics, we can estimate the bandwidth of this 
compensation scheme by relating it to material proper- 
ties of niobium as a type-II superconductor. The field 
suppression factor is determined by how fast the induced 
currents in the superconducting rings respond to changes 
in the external field, which depends on the ring's induc- 
tance (a geometric factor independent of frequency) and 
resistance. Above the first critical field, type-II super- 
conductors exhibit flux pinning, which leads to ac resis- 
tance, but the critical held for niobium is on the order of 
1000 G [l^ll^]. Below the critical field, superconductors 
can still exhibit a frequency-dependent AC resistance as 
described in Ref. [1^. However, for niobium the effect is 
not significant until frequencies up to ~10^^ Hz. There- 
fore at typical bias fields (4 G) and frequencies relevant 
to our qubit (<1 kHz), niobium behaves as a perfect su- 
perconductor and we expect the field suppression factor 
to be the same as that measured under dc. 

Greater reduction can be obtained by optimizing the 
geometry further, for example, by decreasing the distance 
between the plates to 4 mm, but is not implemented 
because of physical constraints in the apparatus. This 
method stabilizes the magnetic field only along the axis 
of the superconducting rings, but since the 4-G bias field 
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FIG. 2: (a) Two superconducting disks, one below and one 
above the trapped ion, stabilize the magnetic field in the z 
direction. (Not to scale.) (b) Magnetic field fluctuation sup- 
pression due to the top disk only (x), the bottom disk only 
(-I-), and both disks (•). When both disks are used, field 
changes are suppressed 50-fold. 



defining the quantization axis is applied in the same di- 
rection, field noise in the x or y direction contributes only 
quadratically to the change in the total field (llj . 



III. MOTIONAL STATE COHERENCE 

The Cirac-Zoller CNOT gate employs superpositions 
of ion motional states as intermediate states during the 
gate, and thus is sensitive to motional decoherence. In 
particular, a high ion heating rate will reduce the gate 
fidelity. An upper bound on the maximum heating rate 
tolerable, /imax, can be given by consideration of the total 
time Tgate required for the pulse sequence implementing 
the CNOT gate, together with a design goal for the gate 
error probability Pgatc desired. Assuming that a single 
quantum of change due to heating will cause a gate error, 
then Umax < T'gatc/Pgate- For Tgate ~ 230 ^s (for our 
experiment), a heating rate of ji-max < 40 quanta/s is 
needed to get Pgatc ~ 0.01. 

We measured the heating rate of the trap at the oper- 
ating secular frequency of 27r x 1.32 MHz. The number 
of motional quanta is measured by probing the blue and 
red sidebands of the S D transition using the shelving 
technique, and comparing the ratio of shelving proba- 
bility on each sideband [29]. The heating rate is deter- 
mined by varying the delay before readout and comparing 
the number of quanta versus delay time. The measured 
heating rate is weakly dependent on the rf voltage and 
dc compensation voltages. Noise on the rf pseudopoten- 
tial can cause heating [30l [3l| , so the ion micromotion is 
minimized using the photon correlation method 32[ . For 
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FIG. 3: (a) Rabi oscillations on the blue sideband. The fitted 
initial contrast is 97.6(3)% and the frequency is 46.7 kHz. (b) 
Ramsey spectroscopy on the blue sideband. The fitted Gaus- 
sian envelope of the decay has time constant T2 — 622(37) /xs. 



correction is done by calculating the shift and account- 
ing for it in the pulse sequencer, following an example in 
NMR [1^. Here, we develop a systematic model of the 
light shifts experienced by a single trapped ion. 

The Stark shift is traditionally a phase shift caused by 
a small change in the transition frequency caused by level 
shifts. In reality, it is a unitary transform involving more 
than just a change of energy levels. We also take this 
into consideration later as a "generalized" Stark shift. 
The model presented here is adapted from well-known 
methods in NMR and included for pedagogical reasons. 
Section llV Al identifles the reference frames useful for dis- 
cussing the single ion in the context of quantum control. 
The generalized Stark shift correction operation is then 
derived as a result of switching between these frames. In 
Sec. lIVBI we apply this method to our single- ion system 
and describe how to calculate the appropriate Stark shift 
correction for any gate in an arbitrary gate sequence. 
Section IIV CI describes the measurement of the ac Stark 
shift and results of the Stark shift correction. 



A. Stark shift on carrier: Simple free-ion model 



more details about the measurements, see Ref. [33[. The 
heating rate can also depend on the trap's proces sing his- 
tory and may vary between temperature cycles [3J|; for 
this trap, the variation is small. In a typical experimental 
run, the rf voltage and dc compensation values are ad- 
justed to minimize the heating rate before the coherence 
time and quantum gate data are taken. Typical heating 
rates obtained in this trap are 4-6 quanta/s, while the 
lowest heating rate measured is 2.1(3) quanta/s. Figure 
^a.) shows Rabi flops on the blue sideband after the ion 
is initialized to the motional ground state with average 
number of quanta n < 0.01. The fitted initial contrast is 
97.6(3)% and the frequency is 46.7 kHz. Motional state 
coherence is demonstrated by performing Ramsey spec- 
troscopy on the blue sideband [FiglJUb)]. The coherence 
time T2* is 622(37) /is. This is comparable to the coher- 
ence time of 660(12) /is of the atomic qubit as measured 
by the same method on the carrier transition. 



IV. STARK SHIFT CORRECTION 

When a two-level atom encoding a qubit is driven off 
resonance, as on a sideband transition, it excites the car- 
rier transition and creates an ac Stark shift. In a real ion 
with multiple levels, additional complication comes from 
other transitions that contribute shifts which are inde- 
pendent of the laser detuning. In the past, correction for 
the Stark shift has been done by using an additional laser 
detuned to the opposite sideband transition to cancel the 
shift 12]. In qubits addressed by a Raman transition, 
this can also be accomplished by changing the power ra- 
tio of the Raman pulses [IHI . In this work, the Stark shift 



There are several useful frames of reference to describe 
the two-level atom model. Consider a single ion at a 
fixed position in free space, interacting with a single- 
mode laser. Let this be described by the laboratory refer- 
ence frame Hamiltonian (with the rotating-wave approx- 
imation) 

Hq = loqcJz + VLcJx cos Lut + flay sinujt (1) 

where ujq is the optical transition frequency, a^, Uy, are 
spin-1/2 operators corresponding to the Pauli matrices 
with eigenvalues ±1/2, and is the Rabi frequency. 

Let the laser be applied at frequency w = + such 
that we may define 



(2) 



as a convenient frame of reference. In the frame of the 
laser, the Hamiltonian is 



(3) 



The frame of reference we wish to use for quantum 
computation (QC frame) is defined by the Hamiltonian 



Thus, if we define a state in this frame as 



(4) 



h{t)}=e+^"<^-'\m) (5) 

where \ip{t)) is the state in the laboratory frame, then we 
find that 



(6) 



assuming that |7(0)) = IV'(O)). 
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The generalized Stark shift correction operation that 
needs to be applied is thus , where 



This is an operator that rotates about an axis 



z + {Q/S)x 



(7) 



(8) 



When the detuning is very large compared with the Rabi 
frequency, the maximum rotation about the x axis, which 
corresponds to a population change, can be bounded by 
fl'^/S^ for a TT rotation about n. For our experimental 
parameters (Sec. IVI Ap . this is less than 1%. Therefore 
the Stark shift is traditionally approximated as a rota- 
tion about the z axis, Rz{d) = e^^"- . We can compute 
what this operation and the rotation angle d would be by 
looking for the Rz closest to R. The angle of rotation of 
the operator ig sjp- ^ fl'^t, while the angle 

of rotation of the operator e^^'^'* is St. Thus, if one ig- 
nored the axes of rotation and treated the first operator 
as if it were also a rotation about z, then the Stark shift 
correction would be a rotation by angle 



(9) 



about z. 



B. Stark shift corrections for arbitrary gate 
sequences 

We now examine a real experimental situation with a 
multilevel ion. To verify our proposed Stark shift correc- 
tion and later to consider the effect of error sources on 
gate fidelity, we simulated gate operations by modeling 
the action of lasers on the full system Hamiltonian in the 
space formed by {|D), \S)} (g) {|0), |1), |2)}. 

Exact simulation of the action of the lasers on the com- 
putational space requires the use of both laser and QC 
frames. The full Hamiltonian is time independent in the 
laser frame, suggesting that gate operations should be 
computed in that frame. The states used in quantum 
computation are defined in the QC frame, where they 
are stationary without the interaction applied. Simula- 
tion of a gate sequence will therefore require frequent 
switching between the frames, for which we define the 
operator C/lqcC*)- 

Computation is performed by moving to the laser 
frame, exponentiating Vl, and moving back to the QC 
frame. For example, a gate performed by application of 
a laser pulse of detuning S, phase 0, starting at time Iq 
for time i, can be computed in the QC frame to be 

U{6,^,t,to) = t/LQc(i + io)e"*''"*/''t/LQc(io)t (10) 

where Vl and ?7lqc depend on the laser detuning, phase, 
and Rabi frequency as well as the trap parameters. Note 



that each laser detuning and trap frequency define a sep- 
arate laser frame. The operator ?7lqc moves between the 
unique QC frame and one of the infinite number of laser 
frames. 



Let 



-20 



be a phase shift on the D states: 



g»0 





(11) 



where /„ is the identity matrix of size n x n. Here the 
2x2 matrix acts on and I3 acts on the mo- 

tional states. Experimentally, this operator is equivalent 
to shifting the laser phase by (/). In a sequence of gates, 
application of such a phase rotation implies shifting the 
laser phases of all subsequent gates. 

In a multilevel atom, there are other transitions that 
are off-resonantly coupled to the laser and contribute to 
additional phase shifts that are detuning independent. In 
our modeling of the Sr+ computation presented here, we 
include the S1/2 O P1/2, S'1/2 ^ P3/2, and D1/2 ^ P3/2 
transitions. The matrix elements of all these transitions, 
which determine the resulting shift, can be calculated as 
in Ref.fH]. 

For gates performed on the carrier transition, since the 
duration of carrier gates is shorter than that of sideband 
gates by the Lamb-Dicke factor 77 (= 0.06 for our case), 
carrier gates take only a small fraction of the total time 
in a typical gate sequence {^2% in the CNOT pulse se- 
quence). Thus we ignore off-resonant coupling to the 
motional sidebands and coupling to far-off-resonant tran- 
sitions. For gates on the sideband transitions, consider 
an interaction with laser detuning S ^ fl, carrier Rabi 
frequency f2, and phase (f) applied for time t starting at 
time to- There are three separate phase shifts that need 
to be cancelled: 

• Stark shift of the ground and excited states can 
be removed by rotating the phase of the |e) state 
hj cj)s = — [5 — V'^^'^i-T?) t, equivalent to applying 
g-*a-j0s^ following the offending gate. Z rotations 
can be performed by changing the phases of all sub- 
sequent laser pulses by (j)s. 

• Resonant excitation of sidebands is applied at a 
frequency Stark-shifted owing to the carrier. The 
laser frame corresponding to that frequency will ro- 
tate with respect to the unshifted states at a rate 
proportional to the Stark shift. To bring the laser 
frame and unshifted states in phase, the laser phase 
has to be shifted by 0/ = (<5 — V S'^ + 51^) ^q. Such 
a phase shift is equivalent to applying e~^'^''^f be- 
fore the gate, and e^'^^'t't after. 

• Off-resonant phase shifts account for approximately 
10% of the total Stark shift in Sr+. Let Aq be 
a constant factor to account for these off-resonant 
phase shifts. 

Define the carrier gate Uc, sideband gate Um, and 
phase correction A as follows. Along with the gate time 
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t and gate starting time to, these variables contain all 
the information relevant to calculating the required Stark 
shift correction. 



a) 



Uc ^ U 

Um - C/0(-A(i + io))C/C/0(Aio) 
A 



5- V^M^ 



(12) 
(13) 

(14) 



From the definition of U and properties of the exponen- 
tial function, it can be shown that the phase correction 
on the nth gate in an arbitrary gate sequence is 



n-l 



J7" = C/0 ( - ^ A'f ) [/ ( (5, </> + A"t^' - ^ AH\ r, 



(15) 

This phase correction consists of the appropriate correc- 
tion for that particular gate plus a global phase, the sum 
of all phase corrections applied to previous gates. In our 
pulse sequencer [23] , the global time and global phase are 
kept as internal registers, and are used to calculate the 
appropriate phase correction every time the qubit laser 
phase is set during a pulse sequence. 



C. Results 

Ramsey spectroscopy on the blue sideband can be used 
to characterize the effectiveness of the Stark shift correc- 
tion. Using the methods described in Ref . . we mea- 
sure the ac Stark shift for various detunings and compare 
to the theoretical model. Figure SJb) shows the typical 
oscillation in the shelving probability P{D) when a pulse 
detuned from the SoD transition (Stark pulse) of vary- 
ing duration is applied. The ac Stark shift is given by 
the oscillation frequency. This shift is measured for sev- 
eral values of detuning and is shown in FigUJ^c) along 
with a one-parameter fit to A/x + b, where the fixed pa- 
rameter is A = fi^/(2ci;sec), the detuning-dependent shift 
to first order, and b is the detuning-independent Stark 
shift caused by farther-off-resonant transitions. The fit- 
ted offset is b — — 27r x 0.5(1) kHz, in agreement with 
—2n X 0.50 kHz predicted by theory. 

The effectiveness of AC Stark shift compensation was 
evaluated by performing Ramsey spectroscopy on the 
sideband and varying the delay time between the two 
pulses, with both Ramsey pulses shifted by 7r/2. In the 
absence of uncorrected Stark shifts, the expected P{D) 
is 1/2 for all delay times. Figure |4l[d) shows the result of 
such a measurement. Here, the secular frequency is deter- 
mined by taking a spectrum and fitting to the sideband; 
then the Ramsey sequence is performed. The parameter 
(j)s is fixed in the experiment control hardware while Ao 
is tuned such that P{D) is maintained near 1/2. From 
the slope of FigHJ^d), we estimate the residual Stark shift 
to be 2tt X 24(20) Hz. 
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FIG. 4: (a) Pulse sequence used to measure ac Stark shift, 
from Ref. . The delay between two carrier 7r/2 pulses is 
Tfixcd = 230/iS while the length r of the Stark pulse is var- 
ied, (b) Typical measurement with the pulse sequence in (a). 
The ac Stark shift is given by the oscillation frequency of the 
shelving probability. Here, the fixed laser detuning is 1.05tJaoc 
which gives an ac Stark shift of 27r x 5.50(4) kHz. (c) Mea- 
sured ac Stark shift as a function of detuning in units of ljscc , 
fitted to the Stark shift model with the detuning-independent 
shift as a free parameter, (d) Ramsey spectroscopy on the 
sideband, demonstrating compensation of Stark shift. 



V. QUANTUM PROCESS TOMOGRAPHY ON 
A SINGLE ION 

With N — 2 qubits encoded in a single ion and meth- 
ods of coupling and controlling these states, a Cirac- 
ZoUer CNOT gate can be implemented [13]. The CNOT 
gate is universal in that all quantum operations can be 
decomposed into single-qubit operations and the CNOT 
gate, and is thus of interest for implementing quantum 
information processing in ion traps. To evaluate the per- 
formance of such a gate, we prepare the system in a set 
of basis states that spans the space of 2^ x 2^ density 
matrices and perform a set of measurements that com- 
pletely specifies the resulting state (state tomography). 
Quantum process tomography (QPT) is performed on 
the two qubits to construct the process matrix, allow- 
ing a full characterization of the gate. Section IV Al gives 
a brief summary of state tomography using conditional 
measurements. Section IV Bl describes a minimal set of 
available measurements and operations in this two-qubit 
single-ion system necessary for QPT. Sections IV CI and 
IV Dl list the pulse sequences for preparing all basis states 
and measuring the outcome. Section|VE] briefly describes 
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the construction of the process matrix that fully charac- 
terizes the gate from these measurements. 



A. Two-qubit state tomography for one ion 

State tomography on the single-ion system of atomic 
and motional qubits requires a nontrivial set of opera- 
tions, since a single qubit rotation on the motional qubit 
cannot be realized directly except by first swapping it 
with the internal state, performing the desired gate, then 
swapping back. The swap operation is complicated since 
the most straightforward physical operations, red- and 
blue-sideband pulses, generally take the system out of 
the computational space, and into higher-order motional 
states such as |2) [3^]. For the CNOT gate, a set of 
composite pulse sequences can keep the system in the 
computational space. But if the goal is measurement 
of the two-qubit state space rather than the realization 
of a coherent operation, an alternative approach can be 
employed. A sequence of measurements, with the second 
conditioned on the results of the first, can suffice to allow 
full state tomography on the two-qubit atomic-|-motional 
state space. This is an extension of the single-ion tomog- 
raphy technique described in [soj . 

The conditional measurement sequence is as follows. 
First we apply an optional tt pulse on the carrier transi- 
tion; then the internal atomic state is measured by flu- 
orescence detection. When this measurement scatters 
photons, it provides information about the internal state 
only and the motional state information is lost. When 
this first measurement does not scatter photons, a tt 
pulse is applied on the blue-sideband transition, which 
allows measurement of the population in the state pairs 
{|S'0),|S'l)} or {|D0),|D1)}, depending on whether the 
initial carrier tt pulse was applied or not. Two measure- 
ments, with and without the carrier pulse, are sufficient 
to determine the population in all four states. 



basis \D0, Dl, D2, SO, SI, S2) are 

^^"(oi)®^3- Po^h-Ps (16) 

The available unitary operations are as follows 

• Rx{d), Ry{ff}- Single qubit (carrier) rotations on 
the { IS*), ID)} qubit. 

• {6) , {6): Blue-sideband rotations, connecting 
{|50),|D1)} and {\Sl),\D2)) (neglecting higher- 
order vibrational modes). 9 is the rotation angle 
on the {\Sff}, \Dl)} manifold. 

• Red-sideband rotations can be defined similarly, 
but are actually not necessary for construction of a 
complete measurement set. 

Explicitly, these rotation matrices are defined as follows: 

Rx{e) = exp[-i0(CT^ (g) /g)] 

Ry{e) ^ cyiY>[-id{(Jy ® h)] 

R+{9) = exp[9{a+(g}a'' - a- (E)a)/2] 

R+{e) = cxp[-i6l(CT+®a^+CT_®a)/2] (17) 

where and a are the creation and annilation operators 
in the Jaynes-Cummings Hamiltonian. 

C. State preparation 

For every measurement sequence, the ion is initialized 
to the state ^I^o = \S0). The sequences of operations 
listed in Table |T] generates the 16 input states that span 
the space of 4x4 density matrices created from the prod- 
uct states \D0,D1,S0,S1). 

D. Complete basis of measurements 



B. Process tomography: Operator definitions 

The state tomographic measurement just described 
measures state populations only, which are the diagonal 
elements of the full density matrix. Relative phases be- 
tween qubit states, which determine coherence properties 
of the state, are also needed in order to perform complete 
process tomography. The phases can be obtained by ap- 
propriate rotations of the qubits prior to measurement. 
Here we define the measurement and rotation operators 
for the sections following. 

The single available measurement is the usual fluores- 
cence detection, which is a projective measurement into 
the \S) state, denoted Ps- Let Pn denote a projection 
into the \D) state. The matrices for Ps and Pd in the 



The following is a procedure for performing complete 
state tomography of the two-qubit {15), \D)} (g) {|0), |1)} 
state of a single ion, using the measurements and opera- 
tions in Sec IV Bl This is a generalization of the method 
used to measure just the diagonal elements of the den- 
sity matrix. There are two kinds of measurement used; 
we call them Mjj and Mjjv- 

Mjj involves performing a unitary operation U on the 
input state and then projecting into the \S) subspace Ps- 
This is described by the measurement operator 

Mu{U) = U^PsU (18) 

Typically, U will be a rotation in the {15*), 1 1?)} subspace, 
implemented by a carrier transition pulse. 

Muv involves first performing a unitary operation U 
on the input state and making a measurement to detect 
fluorescence, which is equivalent to projecting to the \S) 
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Operations applied to \l/o 


State 






\D0) 


*(2) 


R+{-7v) 


\D1) 


*(3) 


I 


ISO) 


*(4) 


Ry{TT)Rt{TT) 


\Sl) 


*(5) 


Rt{Tv)Ry{-n/2) 


{\D0) + \D1))/V2 


*(6) 


R+{-Tv)Ry{-TT/2) 


i\DO) +i\Dl))/V2 


*(7) 


Ry(-n/2) 


{\D0) +\S0))/V2 


*(8) 


7?47r/2) 


{\D0) +i\S0))/V2 


*(9) 


7?^(-7r)i?+(-7r/2) 


{\D0) +\S1))/V2 


^(lO) 


Ry{^iT)R+{n/2) 


{\D0) +i\Sl))/V2 


*(11) 


Rt{Tv/2) 


{\D1) +\S0))/V2 


*(12) 


R+{7v/2) 


{\Dl)+i\S0))/V2 


*(13) 


Ry{7v/2)R+{7T) 


{\D1) + \S1))/V2 


*(14) 


R4-7v/2)Rt{~-n) 


{\D1) +i\Sl))/V2 


*(15) 


Ry{-7v)Rt{-7r)Ry{7T/2) 


{\S0) + \S1))/V2 


*(16) 


Ry{-n)R+{n)Ry{7v/2) 


{\so) + i\siy)/V2 



TABLE I: State preparation operations. 

subspace. Since \D) is long lived, this projection leaves 
the {|I?0), l-Dl), . . .} subspace undisturbed, but motional 
state information is lost if the ion is in state \S,n). If no 
fluorescence is detected, the postmeasurement state is 
PopPo- Conditioned on the first measurement return- 
ing \D) (no fluorescence), a unitary transform V is per- 
formed, and finally another into the 15*) subspace Ps- 
If the first measurement returns fluorescence, the mea- 
surement sequence stops, in which case only information 
about the atomic state is obtained. Mijy is described by 
the measurement operator 

Muv{U, V) = U^PdV'PsVPdU (19) 

Typically, U will be a rotation in the {jS"), jl?)} subspace, 
while V will be one or more rotations on the carrier and 
the red or blue sideband. 

The measurements listed in Table |lT] provide a com- 
plete basis of observables from which the full density ma- 
trix p can be reconstructed, assuming that p is initially 
in only the two-qubit computational subspace. These 
measurement observables are linearly independent. 

The relationship between measurements and the den- 
sity matrix can be expressed by a matrix A with elements 



A, 



(20) 



The full density matrix p can be reconstructed as: 

p = I (21) 



where mj is the result of measurement Mj . 



E. Construction of the process matrix 

A quantum gate including all error sources can be rep- 
resented by the operation £{p), which can be written in 



Ml Muil) 

M2 Muvil, i?+(7r)) 

Ma MuviRyiT^), RtM) 

Mi MuiRyin/2)) 

Ms Mu{R.{-k/2)) 

Me Muv{I, Ry{TT/2)R+{7v/2)) 

Mr MuviRy{7T), Ryin/2)R+iTr/2)) 

Ms Muv{Rvin/2), Ry(n/2)R+ {7t/2)) 

Mg MuviRy{n/2), R4n/2)R+{-K/2)) 

Mio Muvil, R^{-k/2)R+{tx/2)) 

Mil Muv{R^{ti), R^{h/2)R+{h/2)) 

Mi2 Muv{R^{ti/2), R47v/2)R+{7v/2)) 

Mi3 MuviR^-^m, Ry{-K/2)R+{Tx/2)) 

Mi4 Muv{Ry{-n/2), R+{7v/2)) 

Mi5 Muv{R47v/2), R+{7v/2)) 

TABLE II: State measurement functions, 
the operator sum representation as 



(22) 



where p is the input state and Ei is a basis of the set of 
operators on the state space. The process matrix Xmn 
contains the full gate information. For two qubits, the 
state space is spanned by 16 basis states, and 16^ ele- 
ments define the x-matrix, although it only has 16x15 
independent degrees of freedom because of normalization. 
This is reflected in the fact that only 15 measurements 
are needed. The x niatrix can be obtained by inverting 
the above relation. To avoid unphysical results (namely, 
a non-positive-semidefinite p, Tr(p^) > 1) caused by sta- 
tistical quantum error in the experiment, a maximum- 
likelihood estimation algorithm [4X)\ is employed to deter- 
mine the physical operation £ that most likely generated 
the measured data. An alternate, iterative algorithm is 
presented in Ref . [ilj . 



VI. SINGLE-ION CNOT GATE 

The CNOT gate is implemented with the pulse se- 
quence described in Ref. [42| . The optical transition is 
the control qubit, and the motional ground and first ex- 
cited states are used as the target qubit. In the product 
basis {1130), \D1), \S0), 15*1)}, the unitary matrix imple- 
mented is 



f I o\ 

0-100 

1 (23) 
V -10/ 

i(-iy (g)Z + Z(E)I + Z(g)Z + iY(E)Z). 



U 
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This differs from the ideal CNOT matrix by only single- 
qubit phase shifts. Section IVI Al describes the achieved 
gate fidelities and Sec. IVI Bl discusses the major known 
error sources that compromise gate fidelity. 



Gate 


Fr, (%) 


-Fmcan (%) 


Identity 

CNOT 

CNOTx2 


90(1) 
85(1) 
81(1) 


94(3) 
91(5) 
89(6) 



TABLE III: Measured gate fidelities for the identity gate, the 
single CNOT gate, and two concatenated CNOT gates. 



A. Gate performance 



Quantum process tomography was carried out to eval- 
uate the performance of various gates on the two qubits 
of a single ion. The ion in its motional and atomic ground 
state is initialized to one of the 16 input states in Table 
m Then the gate is applied, and the output state is de- 
termined by making all of the measurements listed in Ta- 
ble |lll The longest duration of the full measurement se- 
quence (excluding the gate) is 610 fis, and a single CNOT 
gate takes 230 /zs. These durations are determined by the 
Rabi frequency on the carrier 17 = 27rxl25 kHz and on 
the sideband Hbsb — 27rx7.7 kHz, and the secular fre- 
quency Wscc = 27rx 1.32 MHz. The resulting x matrix for 
the CNOT gate is shown in FiglS] 

We evaluate the performance of the identity gate (all 
preparation and measurement sequences performed with 
no gate in between), the single CNOT gate, and two 
concatenated CNOT gates (CN0Tx2). The results are 
shown in Table Hill The process fidelity is defined as Fp = 
Tr(xidXcxpt), where Xid is the ideal x matrix calculated 
with the ideal unitary operation U, and Xoxpt is experi- 
mentally obtained using maximum-likelihood estimation. 
We also calculate the mean fidelity F^^an, based on the 
overlap between the expected and measured density ma- 
trices, Tr(pidPoxpt), averaged over all prepared and mea- 
sured basis states, as in Ref . [isj . Fp characterizes the 
process matrix whereas Fmean is a more direct measure of 
the gate performance. There exists a simple relation ship 
between the two measures, Fmean = idFp + l)/{d+l) [4J|, 
which is consistent with the independently calculated val- 
ues for our data. Error bars on Fp are calculated from 
quantum projection noise using Monte Carlo methods 
[43]. The large error bars on Fmcan occur because cer- 
tain measured basis states consistently have a higher or 
lower overlap with the ideal states. In general, states that 
involve multiple pulses to create entanglement are more 
susceptible to error and therefore have a lower fidelity 
than states that are closer to pure states. The pulse se- 
quence for some states essentially performs a CNOT gate 
to create and remove entanglement; thus imperfect state 
preparation and measurement contributes significantly to 
the overall infidelity. Using the data for 0, 1, and 2 gates, 
we can estimate the fidelity of a single CNOT gate nor- 
malized with respect to the overall fidelity of the state 
preparation and measurement steps. Assuming that the 
fidelity of the nth gate is F^ = Fi{Fg)^, where Fi is the 
preparation and measurement fidelity, the fitted fidelity 
per gate, Fg, is 95%. 



Error source 


Magnitude 


Approx. contribution 


Off-resonant excitations 


1% 


10% 


Laser frequency fluctuations 


300 Hz 


5% 


Laser intensity fluctuations 


1% 


1% 


Total 




15% 



TABLE IV: Error budget listing the major sources of errors 
on the process fidelity of the single CNOT gate, obtained by 
simulation. Each error source is assumed to be independent. 
The total error is calculated as the product of individual er- 
rors. 



B. Error sources 

A number of possible error sources and their contri- 
butions to the process fidelity of the single CNOT gate 
are listed in Table I IVI To estimate and understand er- 
ror sources, we simulated the full system evolution in 
the (2 atomic state) x (3 motional state) manifold using 
the exact Hamiltonian, including Stark shift and tomo- 
graphic measurements. The magnitude of each source is 
measured independently and then added to the simulated 
pulse sequence. Laser frequency fluctuation is assumed 
to be the primary cause of decoherence and is measured 
by observing the decay of Ramsey fringes on the carrier 
transition. The frequency fluctuation is simulated as a 
random variable on the laser frequency which grows in 
amplitude over time, and accounted for via Monte Carlo 
techniques. Laser intensity fluctuations are measured di- 
rectly with a photodiode. On short time scales compara- 
ble to the length of the gate, the fluctuations are ^ 0.1% 
peak to peak; on longer time scales, up to 1% drifts are 
observed. Both of these effects are accounted for in the 
simulation. Off-resonant excitations are automatically 
included in the model of the full Hamiltonian. The effect 
can be removed from the simulation if decoherence is not 
included and the simulated pulses are of arbitrarily long 
lengths, which is equivalent to reducing the laser inten- 
sity. The resulting x matrix and fidelity agree well with 
the measured results, indicating that the observed fidelity 
is well understood in terms of technical limitations. 

Off-resonant excitations, caused by the square pulse 
shape used to address all transitions, is expected to be 
the largest source of error, as previous work has found 
[4^. Square pulses on the blue-sideband transition con- 
tain many higher harmonics, which causes residual ex- 
citation of the carrier transition. The carrier transition 
oscillations caused by this can be measured directly, av- 
eraged over many scans because of their small amplitude. 
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(a) 



Abs(ZcNOT) 



(b) 



(c) 



'"^tecNoi) 






FIG. 5: Process tomography on the CNOT gate, (a), (b), and (c) show the absolute, real, and imaginary parts of the x matrix, 
respectively. 



Although the measured amount of off-resonant excitation 
is small (^1%) for the laser intensity and secular frequen- 
cies used for our gates, both our simulations and previous 
work 46] have found that up to 10% improvement in gate 
fidelity can be gained by implementing amplitude pulse 
shaping. 



VII. CONCLUSION 

In summary, we have developed a cryogenic, micro- 
fabricated ion trap system and demonstrated coherent 
control of a single ion. The cryogenic environment sup- 
presses anomalous heating of the motional state, as well 
as enableing the use of a compact form of magnetic field 
stabilization using superconducting rings. We perform 
Stark shift correction in the pulse sequencer, removing 
the requirement for a separate laser path and acousto- 
optical modulator. A complete set of pulse sequences 
for performing quantum process tomography on a single 
ion's atomic and motional state is implemented. These 
components are sufficient to perform a CNOT gate on 
the atomic and motional state of a single ion. It is ex- 
pected that amplitude pulse shaping would further im- 



prove the gate fidelity. These techniques, realized in a rel- 
atively simple experimental system, make the single ion 
a possible tool for studying other interesting quantum- 
mechanical systems. 

The control techniques and the CNOT gate demon- 
strated in this work focus on a single ion and do not 
constitute a universal gate set for scalable quantum com- 
putation. However, the additional requirements for such 
a universal two-ion gate, including individual addressing 
[47'] and readout of two ions, have been realized in tra- 
ditional 3D Paul traps as well as other surface trap ex- 
periments, and are not expected to pose significant chal- 
lenges. The microfabricated surface-electrode ion trap 
operated in a cryogenic environment thus offers a viable 
option for realizing a large-scale quantum processor. 
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